Time-resolved homodyne characterization of individual quadrature-entangled pulses 



Jerome Wenger, Alexei Ourjoumtsev, Rosa Tualle-Brouri and Philippe Grangier 
Laboratoire Charles Fabry de I'Institut d'Optique, CNRS UMR 8501, F-91403 Orsay, France. 
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entangled pulses. Quantum correlations equivalent to about 2.5 dB squeezing are efficiently and 
easily reached using the non-degenerate parametric amplification of femtosecond pulses through 
a single-pass in a thin (100 fim) potassium niobate crystal. The entangled pulses are then indi- 
vidually sampled to characterize the non-separability and the entropy of formation of the states. 
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entanglement to the time-resolved homodyne detection. This particularity allows for applications 
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I. INTRODUCTION 

Quantum correlations have properties which cannot be 
reproduced by the rules of classical physics 0, Q . More 
specifically, quantum entanglement is now acknowledged 
as a physical resource, which is needed to perform new 
quantum information processing tasks, such as quantum 
teleportation, dense coding, quantum cryptography or 
quantum computation To perform these protocols, 
quantum continuous variables |4( have recently emerged 
as a relevant alternative to discrete-levels systems. This 
is particularly true in the optical domain, where entangle- 
ment of quantum continuous variables provides a tool of 
major importance for developing new quantum commu- 
nication devices, based on homodyne detection of intense 
beams, rather than photon-counting detectors. 

In order to characterize simply quantum communica- 
tion protocols, it is very convenient to consider them as 
the exchange of some symbols, which carry the relevant 
information and can be accessed individually. Hence, an 
important issue for quantum communication is to develop 
communication schemes operating in a pulsed regime, 
that are able to manipulate individually each quantum 
state involved in the exchange. The analysis of such pro- 
tocols is then very easy in terms of information transfers 
d 0|. However, since the landmark experiment by Ou 
et al. 0, El , most of the experiments relying on entan- 
gled quantum continuous variables have been performed 
by using continuous wave (rather than pulsed) entan- 
gled light beams 0, E3, IH IH 13 Q 111 EE 1H d- 
In addition, the characterization of the quantum corre- 
lations was performed in the spectral domain (i.e. by 
using radio- frequency spectrum analysers), rather than 
in the time domain. As a consequence, even if they are 
actually pulsed 0, |2(J , such experiments cannot be eas- 
ily used to implement quantum communication proto- 
cols, because they do not provide access to each individ- 
ual entangled pulse. These difficulties might be solved 
by developing a quantum theory of analog modulation, 
evaluating the information transfer by defining appropri- 
ate time-bandwidth limited modes, and quantizing them. 
But as long as no such theory has been proposed, ad- 



dressing well-separated "symbols" - each associated to 
one quantum state - clearly appears to be much easier, 
both from a practical and a fundamental point of view. 



In this article, we present a new scheme for the gener- 
ation of quadrature-entangled light pulses, which can be 
individually accessed and characterized. This experiment 
is based on the non-degenerate parametric amplification 
of ultrafast (150 fs) pulses through a single pass in a thin 
potassium niobate (KNbOa) crystal. Thanks to the high 
peak power of the femtosecond pulses and the high non- 
linear coefficient of the KNbC>3 crystal, significant quan- 
tum correlations (equivalent to about 2.5 dB squeezing) 
have been reached. The non-separability of the entangled 
pulses is then directly characterized by a time-resolved 
homodyne detection, which samples the quantum prop- 
erties of each individual incoming pulse. 



Recently, our team has developped a similar setup 
which operated in a degenerate configuration to pro- 
duce pulsed squeezed states [2l|. The experiment de- 
scribed here presents some further developments of this 
scheme towards pulsed quadrature entanglement. As said 
above, all our setup operates in the pulsed regime, from 
the source of entanglement to the time-resolved detec- 
tion, so our analysis is completely carried out in the 
time domain, and not in the frequency domain. Since 
the quadrature-entangled pulses are efficiently and eas- 
ily generated through a single pass in a nonlinear crys- 
tal using non-degenerate parametric optical amplification 
(NOPA), the setup provides a simple and compact source 
for pulsed quadrature entanglement. 



This paper is organized as follows : first, we review 
some general results concerning the manipulation and 
the characterization of quadrature-entanglement (section 
ITT}) . Then the experimental setup is described in details 
(section llllf) before presenting the results of our pulsed 
homodyne measurements (section llVfl . 
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II. QUADRATURE ENTANGLEMENT : 
PRINCIPLES AND CHARACTERIZATION 

A. Producing entanglement : the non-degenerate 
parametric amplifier 

A well-known way to generate a two-mode state whose 
quadrature components are entangled is to use an optical 
parametric amplifier in a non-degenerate configuration 
|22l |23| . Let us denote by (A", P) the quadrature com- 
ponents of one mode of the light field. These quantum 
operators follow the commutation relation [X, P] = 2iNo 
and thus the Heisenberg uncertainty A 2 A A 2 P > N 2 , 
where A 2 (.) denotes the variance of the observable (Nq 
is a scaling constant, which equals to the standard shot 
noise variance). The quadrature components (Xa,Pa) 
and (Xb,Pb) of the output modes ("signal" and "idler") 
of a perfect non-degenerate parametric amplifier are 
given by [H|23| : 



(1) 



Xa = coshr XA,in + sinhr Xb, 

Pa = cosh r Pa,™ -sinh rP B ,i, 

Xb = coshr XB t i n + sinhrXx, 

P B = cosh r P B ,in — sinh tPam 



where (X A ,i n , PA,in) an d (X B ,in , Pb ,in) stand for the in- 
put modes and r is the squeezing factor, which depends 
on the strength of the nonlinear interaction. In the fol- 
lowing, we will consider the case where the input modes 
are in the vacuum state. The state generated by NOPA is 
then generally called two-mode squeezed state, as it bears 
some quantum correlations between the quadratures : 



(X A X B ) = -(PaPb) = sinh2riVo 



((X A - X B ) 2 } = ((P A + P B ) 2 ) = 2 , 



'No 



(2) 
(3) 



In the case of infinite nonlinear effects, r — > oo, the quan- 
tum correlations between the quadratures would become 
perfect (e.g. ((Xa — Xb) 2 ) — > 0). This case corre- 
sponds to the famous Einstein, Podolsky and Rosen 1935 
Gedankenexperiment 0. 



B. Characterizing entanglement 
non-separability criterion 



the 



For a given state composed of two modes A and B, a 
relevant question is to determine whether this state car- 
ries some entanglement or not. For quantum continuous 
variables, Duan et al |24| and Simon |25| have indepen- 
dently formulated a sufficient condition for a state to be 
non-separable : 



i DS = i [a 2 (x a 



X b ) + A 2 (P a + Pb)] <2N (4) 



Dealing with a symmetrical Gaussian two-mode state, 
such as the one generated by NOPA, this condition turns 



out to be a necessary and sufficient condition of non- 
separability. Another criterion used to quantify the quan- 
tum correlations between the entangled modes is the 
"Reid-EPR" criterion [22] , which applies to the products 
of the conditional variances : Vx B \x A Vp b \Pa < Nq 2 . It 
is also known that quantum states that are non-classical 
according to the Reid-EPR condition given above will 
automatically be non-separable according to the Duan- 
Simon condition J3J, though the reverse is not true. 
Thus, in the present work, we will use the Duan-Simon 
criterion as the condition for non-separability. 

In the case of the state generated by perfect NOPA, one 
gets Ids = 2e~ 2r Nq. Thus the violation of the Duan- 
Simon condition for separability occurs as soon as r > 
and increases with the strength r of the nonlinear inter- 
action. This suggests to use Ids to quantify the entan- 
glement. However, a relevant measure of entanglement 
has to satisfy some constraints given in [22, 28] . Unfortu- 
nately, Ids does not fulfil these constraints and thus does 
not provide a "true" measure of entanglement. Neverthe- 
less, Ids is a direct witness of non-separability and can 
be used to compute a relevant measure of entanglement 
for Gaussian states, as we will see now. 



Quantifying entanglement : the entropy of 
formation 



For a bipartite pure state, the amount of entanglement 
is uniquely defined by the Von Neumann entropy of one 
of the sub-state [28J. However, for a mixed state, the def- 
inition of the quantity of entanglement is not unique, and 
has led to different theoretical proposals. Among them, 
let us pay some attention to the entropy of formation, in- 
troduced by Wooters in [2j| . This quantity represents the 
number of qubit pairs necessary to prepare the observed 
correlations. An interesting point of this quantity is that 
it has been explicitly calculated in the case of Gaussian 
states by Giedke and coworkers [30|. For a symmetric 
two-mode Gaussian state whose covariance matrix reads 
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the entropy of formation is then given by the formula : 



(5) 



= f ( \J (V — K X )(V — K p ) J 



(G) 



with f(x) = c+(x) log 2 c+(a;) — C-(x) log 2 c_(a;) and 
c ±( x ) — [x^ 1 / 2 ± x 1 / 2 ] 2 /4. Consequently, for a state 
generated by NOPA, the entropy of formation is directly 
linked to the Duan-Simon quantity Ids as noticed in |18| : 



Ey, 



f 



Ids 

2iVn 



(7) 



Phase 9 




FIG. 1: Simplified scheme of the experimental setup used 
to characterize the quadrature entanglement produced during 
single-pass non-degenerate parametric amplification (NOPA). 
The beam (X+, P+) is then measured by a time- resolved ho- 
modyne detection. 



D. Simple experimental measurements of 
quadrature entanglement 




FIG. 2: Experimental setup. SHG second harmonic gener- 
ation; NOPA non-degenerate parametric amplification; PZT 
piezoelectric transducer; PBS polarizing beamsplitter. 



The method experimentally developed here to charac- 
terize the Gaussian quadrature-entangled states is largely 
inspired from these theoretical results and aims at di- 
rectly measuring the Duan-Simon quantity Ids- The 
principle of our scheme is depicted on Fig. ^ The two en- 
tangled modes generated by NOPA are mixed on a beam- 
splitter of R = 50% reflectivity, while one output mode 
of the beamsplitter is sent to a time-resolved homodyne 
detection for pulse-sensitive quadrature measurements. 

Let 9 denote the relative phase between the entangled 
pulses. When the two entangled modes A,B are mixed 
in phase on the beamsplitter, 9 = 0, the output mode of 
the beamsplitter (AT+,P+) reads : 

1 e r 

X+ = — j= {Xa + Xb) = {XA,in + XB,in) (8) 

Thus the state at the output of the beamplitter is 
a squeezed state, and the corresponding variances in 
quadratures are either above or below the shot noise 
level : A 2 X + = e 2r N Q , A 2 P+ = e~ 2r N . Measuring 
A 2 P + = A 2 (Pa + P b )/2 with the homodyne detection 
gives directly one term to be used in the Duan-Simon 
quantity Ids @ ■ 

If a relative phase of 9 = tt is set between the two en- 
tangled modes A,B, the output mode of the beamsplitter 
(X + ,P + ) then reads : 

1 e~ r 

X+ = ~^{Xa — Xb) = -j={XA,in — XB,in) (9) 



p + = -L(p a -p b) = ^(p a , 



Pf 



This is still a squeezed state, but the squeezed quadra- 
ture has now turned to the X quadrature. Measuring 
A 2 X + = A 2 (X A - X B )/2 with the homodyne detec- 
tion, one gets directly the other term of the Duan-Simon 
quantity Ids- Thus by mixing the entangled pulses on 



a beamsplitter and controlling the relative phase 9 be- 
tween them and between the local oscillator used in the 
homodyne detection, one can measure all the relevant pa- 
rameters to express the Duan-Simon quantity Ids- More 
generally, one can easily show that for any phase 9 the 
quadrature sin(|) X + + cos(|) P + is squeezed while the 
quadrature cos(|) X + — sin(|) P + is anti-squeezed. Thus 
dephasing the EPR beams by 9 before the recombina- 
tion makes the phase-space uncertainty ellipse turn by 
an angle of 9/2, but it keeps the ellipticity constant. 

One interesting feature of this characterization method 
is that it provides simply the necessary quantity Ids to 
check for the non-separability from Eq.Q) and to quan- 
tify the entanglement from Eq.JJJ. Only one homodyne 
detection setup is needed, which greatly simplifies the 
experiment. This method appears also quite similar to 
the setups used for dense coding experiments [TfT . Isil l32l | . 
Let us now turn to the physical implementation of these 
principles in the pulsed regime. 



III. EXPERIMENTAL SET-UP 

We have recently developed a new scheme for pulsed 
squeezed light generation |2l|. This experiment is based 
on the degenerate parametric amplification of femtosec- 
ond pulses. Here we extend this experiment and pass to 
the non-degenerate regime of parametric amplification. 

The experimental scheme is presented on Fig. [5] 
The initial pulses are obtained from a cavity-dumped 
titanium-sapphire laser (Tiger-CD, Time-Bandwidth 
Products), delivering nearly Fourier-transform limited 
pulses at 846 nm, with a duration of 150 fs, an en- 
ergy of 40 nJ, and a pulse repetition rate of 780 kHz. 
These pulses are frequency doubled in a single pass 
through a thin (100 /im) crystal of potassium niobate 
(KNbOs). This crystal is set inside a small vacuum 
chamber and peltier-cooled down to about — 14°C to 
obtain non-critical (90 degrees) type-I phase-matching 
for second harmonic generation (SHG) at 846 nm. The 
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choice for a thin crystal length allows for a wide phase- 
matching bandwidth and avoids the conditions of large 
group-velocity mismatch, contrary to the previously re- 
ported use of thick KNb0 3 crystals [H El- Even for 
the short interaction length used here, KNbC>3 proved 
to be suitable to our applications thanks to its high 
non linear coefficient (about 12 pm/V) and non-critical 
phase-matching. Typically, the SHG efficiency obtained 
is about 28% (corrected from losses). 

A small fraction (1%) of the fundamental beam is 
taken out to serve as a probe to study classical para- 
metric amplification occurring in a similar KNbC>3 crys- 
tal used in a single-pass type-I spectrally degenerate but 
spatially non-degenerate configuration. This spatial non- 
degeneracy is obtained by shifting horizontally the probe 
beam before focussing inside the crystal. Thanks to the 
focussing lens, this shift is transformed into an angu- 
lar difference between the propagation directions of the 
pump and the probe beam within the crystal (see Fig.[3J). 
For small angular differences of about 3° within the crys- 
tal, the phase-matching condition can still be fulfilled by 
lowering the temperature of the crystal by a few degrees. 
The average power of the probe beam after the nonlinear 
interaction is directly measured by a standard Silicium 
photodiode, which allows to estimate for the parametric 
gain. For a well-chosen angular shift so that the probe 
and the pump wave vectors do not overlap, the paramet- 
ric gain does not depend of the relative phase between 
the pump and the probe beam. This condition indicates 
that the non-degenerate configuration is strictly reached 
other the whole spatial extend of the beams. We have 
experimentally optimised the overlap between the pump 
and the probe beam to maximise the amplification gain. 
For the best configuration, the probe waist was set to be 
about y/2 times smaller than the pump waist inside the 
crystal. Then, the best amplification obtained was 1.24. 
This gain was checked to be independent of the probe 
phase and average power. 

The two beams ("signal" and "idler") resulting from 
the classical parametric amplification of the probe are 
mixed on a beamsplitter of R = 50% reflectivity, in order 
to implement the entanglement characterization method 
described in section Hi Dl One output mode of the beam- 
splitter is then sent to the balanced homodyne detection, 
and interferes with the local oscillator beam (LO). The 
homodyne detection is set to be directly sensitive to the 
incoming pulse distribution in the time domain. For each 
pulse, the fast acquisition board (National Instruments 
PCI-6111E) samples one value of th e sig nal quadrature 
in phase with the local oscillator [y, 1.351 l3q . The data 
presented below are obtained directly from these indi- 
vidual pulse measurements. Pulsed homodyning is tech- 
nically more challenging than frequency-resolved homo- 
dyning because low-frequency noises cannot be filtered 
out. Each arm of the detection has to be carefully bal- 
anced (with a typical rejection better than 10 -4 ) even 
for ultra-low frequency noises. By blocking the squeezed 
beam, the detection was checked to be shot-noise limited 
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FIG. 3: Recorded noise pulses while linearly scanning the LO 
phase, when the entangled beams are mixed in phase (9 — 0) 
at the 50 — 50 beamsplitter. Each dot corresponds to the 
measurement of one incoming pulse. 



in the time domain, showing a linear dependence between 
LO power and the noise variance up to 2.5 10 8 photons 
per pulse at a repetition rate of 780 kHz and in the fem- 
tosecond regime. The electronic noise was low enough to 
ensure a ratio larger than 11 dB between shot noise and 
electronic noise variances. 



IV. PULSED HOMODYNE MEASUREMENTS 

A. Characterization of the non-separability 
criterion 

The probe beam being blocked, the amplifier operates 
in a pulsed parametric down-conversion regime to gen- 
erate a two-mode squeezed state. A first piezo-electric 
transducer allows to finely control the relative phase 9 
between the two entangled pulses before the recombina- 
tion beamsplitter. A second piezo transducer is also used 
to control the relative phase of the local oscillator, which 
commands the quadrature component measured by the 
homodyne detection. Fig. [3] displays the recorded noise 
pulses while scanning the local oscillator phase, when the 
quadrature-entangled beams are mixed in phase, 9 = 0. 
Fig. El presents the corresponding quadrature variance 
while scanning the local oscillator phase, when the en- 
tangled beams are in phase, 9 = (a), or dephased by 
6 = 7r (b) at the 50 — 50 beamsplitter. As expected for 
the squeezed states generated by recombining the entan- 
gled pulses, the measured noise variance passes below the 
shot noise level (SNL) for some phase values of the lo- 
cal oscillator. The measured correlation variance (with 
no correction) is 0.70 Nq (-1.55 dB below the shot noise 
level SNL), while the corresponding anti-squeezed vari- 
ance is 1.96 N (2.92 dB above SNL). 

To characterize the entanglement of the two-mode 
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FIG. 4: Quadrature noise variance of the recombined beam 
(X+,P+) (plotted in a linear scale and computed over blocks 
of 2,500 samples) while linearly scanning the LO phase, to- 
gether with the shot noise level (SNL) and the electronic noise 
level. The curve (a) corresponds to the case when the en- 
tangled beams are mixed in phase (9 — 0), while they are 
dephased by 9 — ir for the curve (b). 



FIG. 5: Quadrature noise variance (dark curve) of one entan- 
gled beam (plotted in a linear scale and computed over blocks 
of 2,500 samples) while linearly scanning the LO phase. The 
other entangled beam is blocked before the 50 — 50 beamsplit- 
ter. The second light gray curve above the SNL corresponds 
to the reverse situation by blocking the other beam. 



B. Characterization of the covariance matrix 



state produced by NOPA, it appears logical to compen- 
sate for the homodyne detection efficiency. The proce- 
dure to measure this detection efficiency is well estab- 
lished from squeezing experiments |2l|. The overall de- 
tection efficiency is given as r\ = ijtVhVd = 68%, where 
we have independantly measured the overall transmis- 
sion r]T = 93%, the mode-matching visibility r] B = 88% 
(obtained from the interference fringes between the lo- 
cal oscillator and a seed beam set for maximum classical 
amplification), and the detectors efficiency (Hamamatsu 
S3883) t\d = 94.5%. Given this evaluation of the homo- 
dyne efficiency, one can evaluate the correlation variance 
of the pulses before the homodyne detection, which is 
found to be 0.56 N (-2.52 dB below SNL). 

This experiment was repeated several times for differ- 
ent relative phases 9 between the entangled pulses be- 
fore the recombination. This allowed us to check the 
symmetry of the state produced : we measured the 
same correlation variance for 9 = 0, 9 = tt and sev- 
eral other phases (within reasonable statistical errors less 
than 0.01 No, see FigQJ. Consequently, we can state that 
the correlations between the quadratures are equal to 
A 2 (X A -X B )/2 = A 2 (P A +P B )/2 = 0.56 N . TheDuan- 
Simon quantity amounts then to Zds = 1.12 Nq, which is 
clearly below the threshold for separability (2 Nq) given 
by Eq. This results attests for the quantum entan- 
glement (non-separability) of the pulsed state generated 
by NOPA. One may also quantify the entanglement avail- 
able thanks to the entropy of formation given by Eq. JJJ, 
which provides E-p = 0.44 ebit. 



Theoretically, a general Gaussian state is fully charac- 
terized by its mean values of quadratures together with 
its covariance matrix 7, which comprises the second mo- 
ments of the conjugate quadratures X, P. For a two- 
mode state, the general covariance matrix 7 contains 16 
terms. In our case, the two-mode squeezed state emerg- 
ing from NOPA is generated and handled in a symmetric 
way. This point was experimentally checked as discussed 
in the previous subsection. Consequently, the covariance 
matrix for our state may be reduced to the form given 
by Eq. ©, with V N = A 2 X A = A 2 P A = A 2 X B = 
A 2 P B and K X N = K p N = \{X A X B + X B X A ) = 
— t;(P a P b + P B P A ). This means that for the optimal 
choice of quadratures there is no cross-quadrature cor- 
relations. Let us also point out that for Gaussian states 
the covariance matrix can alwa ys b e reduced to this form 
using local unitary operations |24j . 

By blocking one entangled beam before the 50 — 50 
beamsplitter, we have measured the diagonal terms of 
the covariance matrix (up to the transmission of 50% of 
the beamsplitter and the homodyne detection efficiency 
ff). Our results using the time-resolved homodyne de- 
tection are displayed on Fig. [5] For the two entangled 
beams, we have measured a variance of 1.17 ± 0.01 No, 
which can be related to V = 1.50 before the 50 — 50 
beamsplitter (to get this value, we have corrected for the 
50% transmission of the beamsplitter and the homodyne 
detection efficiency 77 = 68%). The off-diagonal terms 
in the covariance matrix (J5J) can be evaluated by notic- 
ing that thanks to our squeezing measurement A 2 X + we 
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get : 



CONCLUSION 



A 2 A% 



A 2 (X A -X B ) 



(10) 



(A 2 X A + A 2 X B - (X A X B + X B X A )) 



With V 



K„ 



: (V-K X )N Q 

1.50 and A 2 X + /N = 0.56, we get directly 
= 0.94. Thus the reconstructed covariance 



matrix of the experimental entangled pulses reads 

/ 1.50 (0) 0.94 (0) ^ 

(0) 1.50 (0) -0.94 

0.94 (0) 1.50 (0) 

V (0) -0.94 (0) 1.50 / 



(11) 



The values set to zero due to symmetry considerations 
are indicated by parenthesis. A similar procedure has 
already been followed by Bowen et al > but in the case 
of continuous- wave light detected around 3.5 or 6.5 MHz 
using a spectrum analyser. 



An efficient setup to generate and characterize fem- 
tosecond quadrature-entangled pulses is presented, using 
a complete time domain analysis. Pulsed quadrature- 
entangled beams are efficiently generated using non- 
degenerate parametric amplification in a single-pass con- 
figuration. This setup leads to correlations equivalent to 
about 2.5 dB squeezing, that have been recorded using a 
time-resolved homodyne detection. The non-separability 
of the states has then been characterized using the Duan- 
Simon criterion with Ids = 1.12 ATo (< 2Ao), and the 
available entanglement has been quantified using the en- 
tropy of formation to E-p = 0.44 ebit. The association 
of this simple and compact source for quadrature entan- 
glement together with the time-resolved homodyne de- 
tection provides all the ressources for future quantum 
communication protocols using continuous variable en- 
tanglement. 
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